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ABSTRACT 



We present a new, simple method for analysing gapped time-series which is 
particularly suited to probe the X-ray power spectra of active galactic nuclei 
(AGN). The method is applied to a four-day observation of the Seyfert 1 
galaxy MCG —6—30—15 with ASCA (Advanced Satellite for Cosmology and 
Astrophysics). MCG —6—30—15 is well-known for rapid, large amplitude X-ray 
variability and was one of the first AGN to show evidence of relativistic effects 
in its Fe-K line profile. In this source, our method probes the power spectrum in 
the range 10 -3 Hz to 5 x 10~ 2 Hz, thus extending the high-frequency coverage 
by a factor of over 50 compared to previous studies. The ASCA data rule out a 
cessation of variability up to 1.6 x 10~ 2 Hz at greater than 90% confidence. The 
power-law index of the power spectrum in this range is consistent with ~ —1.4. 
At 2.1 x 10 -2 Hz we can place an upper limit of 7% on the intrinsic fluctuations 
of the source luminosity at greater than 99% confidence. Above this frequency 
Poisson noise becomes dominant but we can place an upper limit of 10% on 
intrinsic fluctuations up to 5 x 10~ 2 Hz at 98.9% confidence. The data argue 
against orbital modulation as the origin of the variability since it would imply 
upper limits on the mass of the putative black hole of less than 5 x 1O 5 M , and 
therefore super-Eddington luminosities. 



Subject headings: methods: data analysis - galaxies: active - galaxies: individual: 
MCG -6-30-15 - X-rays: galaxies 
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1. Introduction 

The search for characteristic timescales in the X-ray variability of active galactic nuclei 
(AGN) has been a long-standing problem. The X-ray power spectra have been difficult to 
measure due to a lack of the combination of instrument sensitivity and long, uninterrupted 
exposures. So far, EXOSAT has provided the best measurements (see Lawrence & 
Papadakis 1993; Green, McHardy and Lehto 1993; Czerny and Lehto 1997, and references 
therein). Green et al. (1993), from an analysis of the entire EXOSAT AGN database found 
only nine sources with data of sufficient quality to measure the power spectrum, six of these 
in the range 1(T 4 - 5 - 1(T 3 - 5 Hz and the remaining three (NGC 4051, NGC 5548 and MCG 
—6—30—15) in the range 10~ 5 — 10~ 3 Hz. Green et al. (1993) and Lawrence and Papadakis 
(1993) found the power spectra to be featureless with mean power-law indices —1.72 ± 0.52 
and —1.55 ± 0.09 respectively. No breaks in the power spectra were found in the measured 
frequency ranges. However there have been reports of peaks in the EXOSAT X-ray power 
spectra of three AGN; at 1.7 x 10~ 4 Hz and/or 7.9 x 10~ 5 Hz in NGC 4151 (Fiore, Massaro, 
and Perola, 1989); at ~ 2 x 10~ 3 Hz in NGC 5548 (Papadakis and Lawrence, 1993), 
and at ~ 4 x 10~ 4 Hz in NGC 4051 ( Papadakis and Lawrence, 1995; see also Bao and 
Ostgaard 1994). All of these have yet to be confirmed with other missions. The problem 
is that after the demise of EXOSAT and its highly eccentric orbit, the orbits of subsequent 
satellites have not been able to produce such long, essentially uninterrupted data streams. 
Although sophisticated methods have been devised to analyse time-series with gaps (e.g. 
see Merrifield and McHardy, 1994 and references therein), interpretation of the resulting 
power-spectra, considering the complexities introduced by the 'window-function', is not 
st r aight-f orwar d . 

The only other reports of 'characteristic timescales' have come from a long-term 
variability analysis of NGC 4151 which showed evidence, albeit model-dependent, of three 
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timescales, at ~ 5, 300 — 53, 000 s, 14 days and 183 days (see Papadakis and McHardy, 
1995). So far, a high-frequency cut-off in the X-ray power spectra of AGN has not been 
observed. A search for this cut-off is important because it places hard physical constraints 
on the X-ray emission mechanism and/or the size of the region, and ultimately the mass 
of the putative black hole. So far, the best estimate of the minimum X-ray variability 
timescale has come, for only a few bright AGN, from the points where their EXOSAT 
power spectra hit the Poisson noise (Green et al. 1993). These timescales were all greater 
than 10 3 s. In this paper we extend the search for a high-frequency cut-off in the power 
spectrum down to timescales of 20 s (5 x 10~ 2 Hz) in the Seyfert 1 galaxy MCG —6—30—15, 
utilizing high-sensitivity data taken with ASCA, using a new method for analysing gapped 
time-series. 



We define the 'Excess Pair Fraction' ( hereafter EPF) as a function of the bin size, 
T, of a given lightcurve. For each bin size we count the total number of adjacent pairs of 
non-empty bins and then compute the fraction of those pairs in which the number of counts 
in the two bins differ by, or more than, the sum of the square roots of the counts. This 
fraction is EPF(T). In other words, if we observe C\ counts in the first bin of a pair, for the 



second bin to contribute to EPF(T) it must have C2 counts where \C\ — C2I > \[C~\ + \fC2~ 

This has the solution that C2 must not lie inside the range 1 + C\ ± 2\fC[. The second 
bin is used as the first bin of the next pair tested, and so on. This definition of EPF is 
equivalent to saying that in the Gaussian limit, the counts in each bin of a pair must be 
separated by at least the sum of their standard deviations. 

By definition, the EPF characterizes the source variability on the timescale of the 
bin size and nothing in the definition requires uninterrupted data or inventing fictitious 



2. 



The Excess Pair Fraction 
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data in the gaps. All that is required is that there be enough adjacent pairs of bins. It is 
different to a x 2 test against a constant hypothesis as a function of bin size because the \ 2 
test is global in the sense that it searches for a variability trend over the entire observation 
relative to the mean level. On the other hand, the definition of EPF(T) does not refer to 
a mean, but rather, explicitly measures the bin-to-bin variability. For example, if a source 
showed significant variations on 10 s timescales but no longer term trends, EPF(10) would 
be larger than EPF(1000) because the variability over ~ 1000 s would be 'washed out'. 
The definition of EPF(T) is also manifestly different to the excess variance (e.g. Nandra et 
al. 1997) which probes the integrated power spectrum. 

Figure 1 shows the results of several simulations of the EPF using 10 6 bin-pairs for 
each case. The counts for each bin are drawn from a Poisson distribution in which the 
mean count rate of the source changes from /x to (1 + fi])fi between two adjacent bins in 
a pair, where r] is a uniform random deviate between —1 and +1. Curves of EPF(n) are 
shown for various values of / in the range 0.0-1.0 where /i = HqT. For clarity, only the case 
of a constant source (/ = 0) is shown down to fi < 10 (filled circles). As \i tends to zero, 
EPF(n) for all / tends to zero because of the condition in the definition of EPF(fi) which 
excludes bin-pairs containing one or more empty bins. For a constant source, EPF(fi) is 
almost constant when we are well into the Gaussian regime (// 10). 

Figure 2a shows realistic simulations of EPF(T) for a ~ 4.2 day observation of a 
constant source with a count rate of 1.0 ct/s. The time-intervals for data-accumulation 
were identical to the good-time-intervals (GTI) obtained from a real ASCA observation of 
the Seyfert 1 galaxy MCG —6—30—15 (see Figure 3 and §|3|). However, as is evident from 
the definition of EPF(T), the results do not depend at all on the GTI and EPF(T) is 
consistent with the expectation for a constant source observed with no gaps (solid line, 
taken from the / = simulations in Figure 1). As shown in Figure 1, EPF(T) can be 
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extremely sensitive to variability. The 1.0 ct/s constant source in the above simulation 
was modulated with a sine wave of only 3% amplitude (6% peak-to-peak) and a period 
of 367 s. The resulting EPF(T) is shown in Figure 2b. Not only is the variability easily 
detected, but the periodicity makes characteristic imprints on EPF(T). Namely, EPF(T) 
has minima at times equal to the period and its harmonics. In this example, we would be 
able to determine that 367 s is the fundamental period because no minimum is seen at half 
of this period. 

3. Application to MCG -6-30-15 

The Seyfert 1 galaxy MCG —6—30—15 (z = 0.008) has been extensively studied at 
X-ray and other wavelengths. It is one of the first AGN in which evidence for relativistic 
effects on the Fe-K line profile was found (Tanaka et al. 1995). ASCA observed MCG 
—6—30—15 for over 4 days starting 1994, July 23 (these data have already been discussed 
in other contexts by Iwasawa et al. 1996, and Otani et al. 1996). MCG —6—30—15 is 
highly variable in the X-ray band; the maximum 2-10 keV luminosity was over a factor 
of 7 greater than the minimum luminosity during the four-day observation, in the range 
3.8 - 28.2 x 10 42 ergs s^ 1 (H = 50 km Mpc" 1 s" 1 and q Q = 0). The 0.5-10 keV lightcurve 
obtained from the SIS0 instrument is shown in Figure 3, binned at 512 s (see Tanaka et 
al. 1994 for details of the ASCA instrumentation). The reader is referred to Nandra et al. 
(1997) for details of the ASCA data reduction procedures which are very similar to the ones 
adopted here. 

We constructed lightcurves with a range of time-bin sizes from T = 2 s up to T = 1636 
s. All four instruments aboard ASCA were utilized. Partially-exposed bins were rejected. 
Data in the 0.5-10 keV band were used for the two SIS instruments and in the 0.7-10 keV 
band for the two GIS instruments. However, since the time resolution of SIS data is only 4 s, 
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we did not use SIS data for T < 32 s, allowing sufficient time for several CCD-readouts. The 
time resolution of the GIS data is 2 s or better (depending on bit-rate). EPF(T) was then 
computed for each lightcurve. Note that the lightcurves were not background-subtracted. 
Even when the source has its minimum intensity (see Figure 3), the background constitutes 
less than 11 % of the total counts and is not strongly variable (see Gendreau et at. 1995). 

If we average EPF(T) from different instruments we must model EPF(T) 
independently for the different instruments and then average the models in the same 
way as the data. This is because in general the count rates and GTIs may be different 
for data from the different instruments, and EPF{T) is, for a given T, fundamentally a 
function of the total number of counts per bin (//). Figure 4 shows, using the entire MCG 
—6—30—15 observation, EPF(T) averaged over four instruments for T > 32 s and over two 
instruments otherwise. The error bars are obtained from the square root of the number of 
bin-pairs contributing to EPF(T) divided by the total number of bin-pairs tested. For a 
given mean count rate, /i , we can convert the EPF(fi) models in Figure 1 to EPF(T) by 
interpolation, using fi = fioT. The mean count rates were 1.83, 1.49, 0.96, and 1.14 ct/s for 
SISO, SIS1, GIS2 and GIS3 respectively. In Figure 4, the predicted EPF(T) are shown for a 
source in which the bin-to-bin ratio of mean count rates is independent of T and uniformly 
distributed in the range 1 ± /, for / = (constant source), / = 0.07, and 0.20 (dotted 
curves). It can be seen that the measurement errors cause the EPF{T) to hit the Poisson 
noise below T ~ 50 s. At T = 48s (2.1 x 10 -2 Hz) we can say that the source variability is 
less than 7% (/ < 0.07) at greater than 2.80<r, or greater than 99% confidence. At T = 20 
s, (5 x 10 -2 Hz), approximately the smallest bin size obeying Gaussian statistics, we can 
say that / < 0.10 at 2.53a (98.9 % confidence). At T = 1000 s (10~ 3 Hz) the variability is 
~ 20% (/ ~ 0.20). Note that the low-frequency measurement limit for EPF(T) depends 
on the typical duration of an uninterrupted data train (e.g. satellite orbit) and on the total 
duration of the observation. 
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Detailed fitting of the EPF with models of power spectra will be presented elsewhere. 
What is important here is that the variability amplitude decreases steadily and smoothly 
down to at least 64 s (1.6 x 10~ 2 Hz). There are no sharp features or discontinuities in 
EPF(T) in the range ~ 64 - 1000 s (1.6 x 10~ 2 to 10~ 3 Hz ). Note that if variability ceased 
somewhere in the range T ~ 64 — 1000 s it would have been detected. Yet EPF(T = 64 s) 
is ~ 1.78(7 away from the value for a constant source so such a sharp break can be ruled out 
at > 92% confidence. 

The data are in fact compatible with / oc T r . Figure 4 (solid line) shows a comparison 
of the data with a model with Y = 0.435 (normalization 8.2 x 10~ 3 ). These parameters were 
obtained from an 'eyeball fit' and no attempt was made at formal fitting. Note that T is 
related to the power-law index of the power spectrum, a, by T = — (1 + a). This is because 
EPF(T) is a function of the number of variable bins picked at a given frequency interval 
1/T, and thus is a function of the power in that frequency interval. The power spectrum 
is defined as power per unit frequency so the power-law index for fitting EPF(T) must be 
lowered by unity. Thus the data are consistent with a = —1.435 in the range 5 x 10~ 2 to 
10~ 3 Hz, which is compatible with EXOSAT measurements of MCG —6—30—15 at lower 
frequencies (e.g. a = — 1.48±o;]J by Lawrence and Papadakis 1993 and a = — 1. 36^32 by 
Green et al. 1993). 

If the variability that is observed down to T = 48 s (2.1 x 10~ 2 Hz) is due to modulation 
by matter orbiting the putative black hole, then we can infer an upper limit on the 
black-hole mass by requiring that the emitting matter is not located closer than the radius 
for marginally stable circular orbits (i.e. 3 Schwarzschild radii). This condition means 
that the observed period, P, for Keplerian orbits around a Schwarzschild and a maximally 
rotating Kerr black hole is P > 1.05 x 10 4 M 8 s and P > 4.6 x 10 4 M 8 s respectively (M 8 is 
the mass in units of 1O 8 M ). Then P = 48 s implies M 8 < 4.6 x 10~ 3 for a Schwarzschild 
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metric and M 8 < for a maximally rotating Kerr metric. These masses are rather small 
and imply luminosities in excess of the Eddington limit. This argues against the observed 
variability being due to orbital modulation. Note that the upper limits on the black hole 
masses obtained by Green et al. (1993) are in error by a factor of 10 because the numerical 
factor in their equation 6 is incorrect and should be 10 3 and not 10 4 . Also, they arbitrarily 
assume an orbital radius of 5 Schwarzschild radii. These factors, along with the fact that 
the ASCA data detect variability at shorter timescales than the EXOSAT data, explain 
why our upper limit for the black-hole mass in MCG -6-30-15 (under the assumption of 
orbital modulation around a Schwarzschild hole) is a factor of 30 less than the Green et al. 
(1993) estimate of 1.4 x 1O 7 M . 

Clearly, to infer a source size from the fastest observed variability is too simplistic. 
Instead, we have begun, for the first time, to probe the source sub-structure since the rapid 
variability may be produced by localized flares or hotspots. The lower limit on any break 
frequency of 1.6 x 10~ 2 Hz (64 s) represents an upper limit (1.92 x 10 12 cm) only to the size 
of a region responsible for the shortest observed (and low-amplitude) variability, and not 
necessarily to the size of the source region as a whole. 



4. CONCLUSIONS 

We have presented a new method for analysing gapped time-series which is capable of 
probing the power-spectrum of AGN up to ~ 5 x 10~ 2 Hz for sources with count rates of 
the order ~ 1 ct/s or more. We have applied the method to a four-day ASCA observation 
of the Seyfert 1 galaxy MCG —6—30—15 to search for a high-frequency cut-off in the power 
spectrum. The data are consistent with a power-law power spectrum, with no break in 
the range 10~ 3 Hz to 1.6 x 10~ 2 Hz (at a confidence level of 92%). The power-law index 
of the power spectrum in this range is consistent with ~ —1.4. In the range 1.6 x 10~ 2 
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Hz to 5 x 10~ 2 Hz, the fluctuations in the intrinsic mean count rate of the source can be 
constrained to be less than 7% at greater than 99% confidence. It is likely that the observed 
variability is not due to orbital modulation since black hole masses of less than 5 x 1O 5 M 
and super-Eddington luminosities are implied. 



The authors would like to thank Karen Leighly for some extremely useful discussions 
and an anonymous referee for some very helpful criticism. 
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Figure Captions 



Figure 1 

The 'Excess Pair Fraction', EPF(fx), calculated from simulations using 10 6 bin-pairs. The 
mean number of counts in each bin is + frj), where rj is a uniform random deviate 
between —1 and +1 and / represents a fractional variability parameter. The curves are 
shown for several values of / in the range 0.0-1.0 (a value of 0.0 corresponds to a constant 
source). The error bars are obtained from the square root of the number of bin-pairs 
contributing to EPF(/i), divided by the total number of bin-pairs. Note that although the 
error bars are plotted, they are too small to be visible for most of the parameter space. See 
text for details. 

Figure 2 

The 'Excess Pair Fraction', EPF(T), calculated from simulations of a source with mean 
count rate 1.0 ct/s. Good Time Intervals (GTI) were taken from a real, ~ 4.2 day 
ASCA observation of MCG —6—30—15 (see Figure 3). Panel (a) shows the case for a 
constant source and panel (b) shows the same source with a 3% sinusoidal modulation (6% 
peak-to-peak) and a period of 367 s. 

Figure 3 

The 0.5-10 keV SIS0 lightcurve for the extended ASCA observation of MCG -6-30-15. 
The lightcurve is binned at 512 s. The reference time is 23 July 1997 UT 05:56:13. 

Figure 4 

The 'Excess Pair Fraction' (EPF) calculated from the extended ASCA observation of 
MCG -6-30-15 (filled circles). The data from the two SIS instruments (0.5-10 keV) and 
two GIS instruments (0.7-10 keV) were used. For bin sizes less than 32 s only GIS data 
were used, the EPF being averaged for GIS2 and GIS3. Otherwise, the EPF was averaged 
for all four ASCA instruments. The error bars are obtained from the square root of the 
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number of bin-pairs contributing to EPF, divided by the total number of bin-pairs tested. 
The dotted lines show the predicted EPF for f = 0% (constant source), 7% and 20% (as 
indicated), where / is the maximum bin-to-bin fluctuation in the mean count rate of the 
source. The mean number of counts in each bin of width T is /iqT(1 + frj) where i] is a 
uniform random deviate between —1 and +1. The solid curve corresponds to a model with 
/ = 8.2 x io- 3 T+ - 435 . See text for details. 
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